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Abstract. In this paper we study analytically the viscous 'sabra' shell model 
of energy turbulent cascade. We prove the global regularity of solutions and 
show that the shell model has finitely many asymptotic degrees of freedom, 
specifically: a finite dimensional global attractor and globally invariant inertial 
manifolds. Moreover, we establish the existence of exponentially decaying 
energy dissipation range for the sufficiently smooth forcing. 
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1. Introduction 



Shell models of turbulence have attracted interest as useful phenomenological 
models that retain certain features of the Navier-Stokes equations (NSE). Their 
central computational advantage is the parameterization of the fluctuation of a 
turbulent field in each octave of wave numbers A™ < \k n \ < X n+1 by very few repre- 
sentative variables. The octaves of wave numbers are called shells and the variables 
retained are called shell variables. Like in the Fourier representation of NSE, the 
time evolution of the shell variables is governed by an infinite system of coupled 
ordinary differential equations with quadratic nonlincarities, with forcing applied 
to the large scales and viscous dissipation effecting the smaller ones. Because of 
the very reduced number of interactions in each octave of wave numbers, the shell 
models are drastic modification of the original NSE in Fourier space. 

The main object of this work is the "sabra" shell model of turbulence, which was 
introduced in [26]. For other shell models see [3], [12], [13], [31]. A recent review 
of the subject emphasizing the applications of the shell models to the study of the 
energy-cascade mechanism in turbulence can be found in [2] . It is worth noting that 
the results of this article apply equally well to the well-known Gledzer-Okhitani- 
Yamada (GOY) shell model, introduced in [31]. The main difference between shell 
models and the NSE written in Fourier space is the fact that the shell models contain 
only local interaction between the modes. We like to point out that our analysis 
could be performed also in the presence of the long, but finite, range interactions. 

The sabra shell model of turbulence describes the evolution of complex Fourier- 
like components of a scalar velocity field denoted by u n . The associated one- 
dimensional wavenumbers are denoted by fc„, where the discrete index n is referred 
to as the "shell index" . The equations of motion of the sabra shell model of turbu- 
lence have the following form 

= i{ak n+ iu n+2 u* n+l + bk n u n+ iu* n _ 1 - c/j„_iu„_iu„_ 2 ) - vk^Un + fn, (1) 

for n — 1,2,3,..., and the boundary conditions are m_i = u = 0. The wave 
numbers k n are taken to be 



with A > 1 being the shell spacing parameter, and k > 0. Although the equation 
does not capture any geometry, we will consider L — k^ 1 as a fixed typical length 
scale of the model. In an analogy to the Navier-Stokes equations v > represents 
a kinematic viscosity and /„ are the Fourier components of the forcing. 

The three parameters of the model a, b and c are real. In order for the sabra 
shell model to be a system of the hydrodynamic type we require that in the inviscid 
(u = 0) and unforced (/„ = 0, n = 1, 2, 3, . . . ) case the model will have at least one 
quadratic invariant. Requiring conservation of the energy 



leads to the following relation between the parameters of the model, which we will 
refer as an energy conservation condition 



k n = fc A™, 



(2) 



■DC- 




(3) 



a + b + c = 0. 



(4) 
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Moreover, in the inviscid and unforced case the model possesses another quadratic 
invariant 



oo 



w=vr ki 2 . (5) 



n=l 

The physically relevant range of parameters is \a/c\ < 1 (sec [26] for details). 
For — 1 < - < the quantity W is not sign-definite and therefore it is common to 
associate it with the helicity - in an analogy to the 3D turbulence. In that regime 
we can rewrite relation (5) in the form 

oo 

W = £(-l)"fc>„| 2 , (6) 

n=l 

for 

a = log A C - . (7) 

The 2D parameters regime corresponds to < | < 1. In that case the second 
conserved quadratic quantity W is identified with the enstrophy, and we can rewrite 
the expression (5) in the form 

oo 

w = ]Tfc«K| 2 , (8) 

71=1 

where a is also defined by the equation (7). 

The well-known question of global well-posedness of the 3D Navier-Stokes equa- 
tions is a major open problem. In this work we study analytically the shell model 
(1). Specifically, we show global regularity of weak and strong solutions of (1) 
and smooth dependence on the initial data. This is a basic step for studying the 
long-time behavior of this system and for establishing the existence of a finite di- 
mensional global attractor. 

Our next result concerns the rate of the decay of the spectrum (i. e. \u n \) as 
function of the wavenumber k n . Having shown the existence of strong solutions, we 
prove that \u n \ decays as a polynomial in k n for n large enough, depending on the 
physical parameters of the model. However, if we assume that the forcing /„ decays 
exponentially with k n , which is a reasonable assumption, because in most numerical 
simulations the forcing is concentrated in only finitely many shells, then we can show 
more. In fact, following the arguments presented in [21] and [20] we show that in 
this case the solution of (1) belongs to a certain Gevrey class of regularity, namely 
\u n \ also decays exponentially in k n - an evidence of the existence of a dissipation 
range associated with this system, which is an observed feature of turbulent flows 
(see also [11] for similar results concerning 3D turbulence). Similar results on the 
existence of the exponentially decaying dissipation range in the GOY shell models 
of turbulence were obtained previously in [33] by the asymptotic matching between 
the non-linear and viscous terms. Our results set a mathematically rigorous ground 
for that claim. Moreover we show that the exponentially decaying dissipation range 
exists for all values of the parameters of the model, in particular also for c = 0, 
improving the result of [33] . 

Classical theories of turbulence assert that the turbulent flows governed by the 
Navier-Stokes equations have finite number of degrees of freedom (see, e.g., [13], 
[25]). Arguing in the same vein one can state that the sabra shell model with non- 
zero viscosity has finitely many degrees of freedom. Such a physical statement can 
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be interpreted mathematically in more than one way: as finite dimensionality of the 
global attractor, as existence of finitely many determining modes, and as existence 
of finite dimensional inertial manifold. 

Our first results concern the long time behavior of the solutions of the system (1). 
We show that the global attractor of the sabra model has a finite fractal and Haus- 
dorff dimensions. Moreover, we show that the number of the determining modes 
of the sabra shell model equations is also finite. Remarkably, our upper bounds 
for the number of the determining modes and the fractal and Hausdorff dimensions 
of the global attractor of the sabra shell model equation coincide (see Remark f ). 
This is to be contrasted with the present knowledge about Navier-Stokes equations, 
where the two bounds are qualitatively different. In the forthcoming paper we will 
investigate the lower bounds on the Hausdorff dimension of the global attractor. 

Finally, and again in contrast with the present knowledge for the Navier-Stokes 
equations, we show that the sabra shell model equation has an inertial manifold 
M. . Inertial manifolds are finite dimensional Lipschitz globally invariant manifolds 
that attract all bounded sets in the phase space at an exponential rate and, in 
particular, contain the global attractor. We show that M. = graph($>), where $ 
is a C 1 function which slaves the components u n , for all n > N, as a function of 
i u k}k=i' f° r N large enough depending on the physical parameters of the problem, 
i.e. v, /, A, a, b, and c. The reduction of the system (I) to the manifold A4 yields 
a finite dimensional system of ordinary differential equations. This is the ultimate 
and best notion of system reduction that one could hope for. In other words, an 
inertial manifold is an exact rule for parameterizing the large modes (infinite many 
of them) in terms of the low ones (finitely many of them). The concept of inertial 
manifold for nonlinear evolution equations was first introduced in [17] (see also [8], 



2. Preliminaries and Functional Setting 

Following the classical treatment of the NSE, and in order to simplify the nota- 
tion we are going to write the system (1) in the following functional form 



in a Hilbert space H. The linear operator A as well as the bilinear operator B will 
be defined below. In our case, the space H will be the sequences space £ 2 over the 
field of complex numbers C. For every u,v e H, the scalar product (•, •) and the 
corresponding norm | • | are defined as 



We denote by {(t>j}JL 1 the standard canonical orthonormal basis of H, i.e. all the 
entries of 4>j arc zero except at the place j it is equal to I. 

The linear operator A : D(A) — > H is a positive definite, diagonal operator 
defined through its action on the elements of the canonical basis of H by 



[9], [18], [34], [36]). 



— + vAu + B(u, u) = f 
u(0) = u m , 



(9a) 
(9b) 




A(j)j = tffa, 
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where the eigenvalues k 2 satisfy the equation (2). D(A) - the domain of A is a 
dense subset of H. Moreover, it is a Hilbert space, when equipped with the graph 
norm 

\\u\\ D(A) = \Au\, VueD(A). 

Using the fact that A is a positive definite operator, we can define the powers A s 
of A for every sel 

Vw= (ui,tt2,ti3,...), A s u= (k 2s u 1 ,k 2s u 2 ,k' 2s u 3 , . . .). 

Furthermore, we define the spaces 

oo 

V s :=D{A s / 2 ) = {u=(u 1 ,u 2 ,m,...) : ^ kf \uA 2 < oo}, 

j'=i 

which arc Hilbert spaces equipped with the scalar product 

(«, v) a = (A s / 2 u, A s / 2 v), Vu, v e £>(A S / 2 ), 

and the norm |w| 2 = (it, u) s , for every it S Z?(A S / 2 ). Clearly 

V,CV = fiCV.„ Vs>0. 

The case of s = 1 is of a special interest for us. We denote V = D{A 1 / 2 ) a 
Hilbert space equipped with a scalar product 

((u, v)) = {A 1 / 2 u 1 A 1 ' 2 v), Vu, w G £(,4 1/2 ). 

We consider also V' = Z?(A -1 / 2 ) - the dual space of V. We denote by (•, •) the 
action of the functionals from V on elements of V. Note that the following inclusion 
holds 

V C H = H' C V, 

hence the H scalar product of / € H and u E V is the same as the action of / on 
u as a functional in the duality between V' and V 

(/,u) = (/,«), V/eff, VueV. 

Observe also that for every u € D(A) and every d£V one has 

((it, v)) — (Au, v) = (Au, v) . 

Since D(A) is dense in V one can extend the definition of the operator A : V — > V 
in such a way that 

(Au, i>) = ((u,v)), Vit, i> € V. 
In particular, it follows that 

puiiv^iMi, v ue y. (io) 

Before proceeding and defining the bilinear term B, let us introduce the sequence 
analogue of Sobolev functional spaces. 

Definition 1. For 1 < p < oo and m e R we define sequence spaces 

/ OO s 1/p 

u^* :={«=(« 1)U2 ,...)e: ||^ ro/2 «|| P = ( S^CKIH <TO >' 

for 1 < p < oo, and 

w ™<°° : ={u= (ui,U2,...) € : P m/2 «||oo - sup (CKI) < oo}. 

l<n<oo 
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For u £ w m,p we define it's norm 

\\u\\ wm , P = \\A m ' 2 u\\ p , 

where \\-\\ p is the usual norm in the i v sequence space. The special case of p = 2 and 
m > corresponds to the sequence analogue of the classical Sobolev space, which 
we denote by 

h m =w m ' 2 . 

Those spaces are Hilbert with respect to the norm defined above and its corresponding 
inner product. 

The above definition immediately implies that h 1 = V. It is also worth noting 
that V C w 1,oc and the inclusion map is continuous (but not compact) because 

\\u\\ w i,oo = \\A 1 / 2 u\\ 00 <\\A 1 / 2 u\\ 2 = \\u\\. 

The bilinear operator B(u,v) will be defined formally in the following way. Let 
u,v £ H be of the form u = ^2n=i M n^n an< i v — X^raLl v n ( f ) n- 

Then 



OO / 

B(U,V) =~i^2 [ afc n+l«n+2<+l + bk n v n+1 u^_ 1 + 
n=l ^ 



+ ak n -iu n -\v n -2 + bk 

(ii) 

where here again uo — tt_i = vq = V-\ = 0. It is easy to see that our definition of 
B(u,v), together with the energy conservation condition 

a+b+c= 

implies that 

B(u,u) = -iy~] I ak n+1 u n+2 u* n+1 + bk n u n+ iu* n _ x 

n=l ^ ' 

which is consistent with (1). In the sequel we will show that indeed our definition 
of B(u, v) makes sense as an element of H , whenever u £ H and v £ V or u £ V 
and v £ H. For u,v £ H we will also show that B(u, v) makes sense as an element 
of V. 

Proposition 1. 

(i) B : H x V — > H and B : V x H — > H are bounded, bilinear operators. 



Specifically, 
and 

where and 



\B(u,v)\ < CiHMI, Vu£H,v£V, (12) 
|^(u,v)| < C 2 ||«|||vj, Vu£V,v£H, (13) 



C 1 = (H(A- 1 +A) + |6|(A- 1 + 1)), 

C 2 = (2|a|+2A|6|). 
(ii) B : H X H — > V is a bounded bilinear operator and 

\\B{u,v)\\ v > <d\u\\v\, Vu,v£H. (14) 
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(iii) B : H x D{A) — > V is a bounded bilinear operator and for every u G H 
and v e D(A) 

\\B(u,v)\\ <C 3 \u\\M, (15) 

where 

C7 3 = (M(A 3 + A- 3 ) + H(A + A- 2 )). 

(iv) For every u G H , v G V 

Re(B(u,v),v) = 0, (16) 

Proof. To prove the first statement, let u £ H and u G V. i? is a Hilbert space, 
therefore, 

|B(m,u)|= sup \(B(u,v),w)\. 
weH,\w\<i 

Using the form of the eigenvalues k n given in (2) and applying the Cauchy-Schwartz 
inequality we get 



\B(u,v)\ — sup 

weH,\w\<i 



^2 ak n+ iv n+2 u* n+1 w* n + bk n v n+1 w* n u* n _ 1 + 



n=l 



n —\W n V n —\U n —2 



< 



< {\a\{\- 1 + X) + \b\{\~ 1 + 1)) sup MIML^M < Ci|u|||u||, 

w£H,\w\<l 

where the last inequality follows from the fact that ||t;|| „,!,«> < ||t>||. The case of 
u G V and v G H is proved in the similar way. 
As for (ii), let u,v & H and w eV. Then 



\(B(u,v),w}\ 



< 



< 



^2,B n (u,v)w* n 

71=1 
oo 

^ afc n+ iu„ +2 < +1 w;* + Wc„v„ + iu;* w*_ r 



71=1 



+ afc„_iw*u„_iu„_2 + 6fc„_iw*w„_iu n _ 2 



< 



< (laKA- 1 + A) + l&KA- 1 + 1))HU~MM < Ci||HI |«|M. (17) 

Similar calculations can be applied to prove the statement (iii). 

The statement (iv) follows directly from the energy conservation condition (4). 

□ 



Let us fix m G N and consider H m — span{<p\, . . . , <p m } and let P m be an £ 2 
orthogonal projection from £ 2 onto H m . If u G H, then 

P m u = u m = «, 0,0,...), 

where w™ = (u,4> n ). In order to simplify the notation we will write u n instead of 
u" 1 , whenever it will not cause any confusion. We would like to obtain the same 
estimates as above in Proposition 1 for the truncated bilinear term P m B(u m , v m ). 
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Setting Uj = u" 1 = for all j = 0,-1 and j > m we get 



3 

m— 1 



in — j. ^ 

P m B(u m ,v m ) = -i^2 \ ak n+iv n +2u* n+1 (t> n -i + bk n v n+1 u* n _10„+ 

n=2 ^ 

+ afc„_iU„_ii;„_20n+l + ^„_il)„_iU„_2<An+l^ • 

We observe that similar estimates to (14), (18) and (16) hold for the truncated 
bilinear operator P m B(u m ,v m ). Namely, Vu m ,v m e H m 

\\P m B(u m ,v m )\\ v , <Ci|« m ||w TO |. (18) 

and \Ju m ,v m e H 

Re(P m B(u m ,v m ),v m ) = 0. (19) 

3. Existence of Weak Solutions 

In this section we prove the existence of weak solutions to the system (9). Our 
notion of weak solutions is similar to that of the Navier-Stokes equation (see, for 
example, [5], [35]). 

Theorem 2. Let f e L 2 ([0, T), V) and u m e H. There exists 

U eL°°([0,T],!/)nL 2 ([0,T],n (20) 

with 

^i 2 ([0,T],y'), (21) 
satisfying the weak formulation of the equation (9), namely 

^,v^+v (Au, v) + (B(u, u), v) = (/, v) (22a) 
u(jc,0) = u m (x), (22b) 



for every v e V. 

The weak solution u satisfies 

ugC([0,T],H). (23) 

Proof. The proof follows standard techniques for equations of Navier-Stokes type. 
See, for example, ([5], [35] or [36]). We will use the Galerkin procedure to prove 
the global existence and to establish the necessary a priori estimates. 

The Galerkin approximating system of order m for equation (9) is an m - di- 
mensional system of ordinary differential equations 

dii m 

— + vAu m + P rn B{u m , u m ) = P m f (24a) 

u m (0) = P m u m . (24b) 

We observe that the nonlinear term is quadratic in u m , thus by the theory of 
ordinary differential equations, the system (24) possesses a unique solution for a 
short interval of time. Denote by [0,T m ) the maximal interval of existence for 
positive time. We will show later that in fact, T m — T. 
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Let us fix m and the time interval [0, T]. On [0, T m ) we take the inner product 
(-, -)h of the equation (24a) with u m , considering the real part and using (19) we 
get 

~\u m \ 2 + v\\u m \\ 2 = Re(P m f,u m ). 

The right hand side could be estimated using the Cauchy-Schwartz and Young 
inequalities 

\{P m f,u™)\<\\f\\ vl \\u™\\ <^\\f\\l, + V -\\ u ™\\\ 
Plugging it into the equation we get 

±\u™\ 2 + v\\u m \\ 2 < l -\\f\\ 2 v ,. (25) 

Integrating the last equation from to s, where < s < T m , we conclude that 
for every such s, 

\u m (s)\ 2 < \P m u in \ 2 + - f WfWWvdt. 

v Jo 

Thus, limsup s ^ T - |w m (s)| 2 < oo. Therefore, the maximal interval of existence of 
solution of the system (24a) for every m is [0, T) - an interval, where the terms of 
the equation make sense, and in particular, /. If / € L 2 ([0, oo), V), then we can 
take T > arbitrary large and show the existence of weak solutions on [0, T). Now, 
we observe that 

IKI| 2 L~([0,T),i/) < \PmU in \ 2 + l£\\f\\ V ,dt, 

and hence 

The sequence u m lies in the bounded set of L°°([0, T],H). (26) 
Next, integrating (25) over the interval [0, T] we get 

\u m {T)\ 2 + u [ T \\u m (t)\\ 2 dt<\u m (0)\ 2 + - [ T \\f(t)\\ v ,dt, 
Jo v Jo 

and using the fact that |w m (0)| = \P m u m \ < \u m \ we conclude 

The sequence u m lies in the bounded set of L 2 ([0, T], V). (27) 

Finally, in order to apply Aubin's Compactness Lemma (see, e.g., [1], [5]) we 
need to estimate the norm H^jfHIy/- Consider the equation 

dv m 

— = -vAu m - P m B(u m ,u m ) + P m f 

Then, (18) and (26) imply that \\P m B(u m , u m )\\ v , is uniformly bounded for every 
< t < T, hence 

P m B(u m ,u m )eL p ([0,T},V'), 
for every p > 1 and T < oo. Moreover, from the fact that ||-P m /||y/ < ||/||y and 
the condition on / it follows that P m f e L 2 ([0,T], V). 

The relations (10) and (27) imply that the sequence {Au" 1 } is uniformly bounded 
in the space L 2 ([0,T],V). Therefore, using the fact that / e L 2 ([0, T], V) we 
conclude, 

du m 

The sequence lies in the bounded set of L 2 ([0, T], V'). (28) 
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At last, we are able to apply Aubin's Compactness Lemma. According to (26), 
(27) and (28), there exists a subsequence u m , and u £ H satisfying 

u m ' — >u weakly in L 2 ([0,T},V) (29a) 

u m ' — >u strongly in L 2 ([0,T},H) (29b) 

We are left to show that the solution u that we found, satisfies the weak for- 
mulation of the equation (9) together with the initial condition. Let v £ V be an 
arbitrary function. Take the scalar product of (24a) with v and integrate for some 
< t < t < T to get 

(U m '(t),v) - (U m '(t ),v)+V [\{u m '(T),v))dT+ 

+ [ {P m ,B{u m \u m '),v)dT= [ (P m ,f,v)dT. (30) 

The conclusion (29) implies that there exists a set E C [0, T] of Lebesgue measure 
and a subsequence u m of u m such that for every t £ [0, T] \ E, u m (t) converges 
to u{t) weakly in V and strongly in H . As for the third term 

lim f\(u m '(r),v))dT^ [ ((u(r),v))dr, 



I to J tQ 

because of (29a). The same holds for the forcing term. We are left to evaluate the 
nonlinear part. It is enough to check that 
rt 

(P m ,B(u m \u m '),v)- (B(u,u),v)\dT — >0, asm'^oo. (31) 

Let us define Q = Q m i = I — P m < . Then, by definition of B(u, u) and arguments 
similar to (14) we have 



/' 

Jt„ 



\{P rn ,B{u m ' ,u m '),v)- B{u,u) lV )\ = 

OO s 

^2 ( ak n+ iv n+2 ul +1 v* + bk n v n+1 v^ l u* n _ 1 + ak n -iv*u n -iv n - 2 + 

n.—m ^ 



+ bk n -iv*v n -\u n -2 



<C||O w |||Q u | 2 <C||i;||| U -P m m |2 



We see that (31) holds, because of (29b). Passing to the limit in (30) wc finally 
conclude that u is weakly continuous in H and that it satisfies the weak form of 
the equation, namely 

{u{t),v)-{u{t ),v) + y f {{u{T),v))dT+ 

Jto 



+ f {B(u, U ),v)dT = f (f,v)dT. 
Jt Jto 



In order to finish the proof we need to show that u satisfies (23). First, recall 
that u e L 2 ([0,T],VO and ^ e L 2 ([0, T], V). Therefore, it follows that \u{t)\ 2 is 
continuous. Moreover, u £ C w ([0, T], H) - the space of weakly continuous functions 
with values in H. Now, let us fix t £ [0,T] and consider the sequence {t n } C [0,T] 
converging to t. We know that |w(i„)| converges to \u(t)\ and u(t n ) converges 
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weakly to u(t). Then, it follows that u{t n ) converges strongly to u(t), and we 
finally conclude that 

UGC([0,T],H). 

□ 



4. Uniqueness of Weak Solutions 

Next we show that the weak solutions depend continuously on the initial data 
and in particular the solutions are unique. 

Theorem 3. Let u(t), v(t) be two different solutions to the equation (9) on the time 
interval [0, T] with the corresponding initial conditions u m and v m in H. Then, for 
every t G [0, T] we have 

\u(t) -v(t)\ < e K \u m -v in \, (32) 

where 

K = Ci f \\u(s)\\ds. (33) 
Jo 

Proof. Consider u, v - two solutions of the equation (9a) satisfying the initial con- 
ditions u(0) = u m and v(0) = v m . The difference w = u — v satisfies the equation 

^ + vAw + B(u, w) + B{w, u) - B(w, w) = 0, 

with the initial condition w(0) = w m = u m — v m . 

According to Theorem 2, ^ e L 2 ([0,T),V) and w e L 2 ([Q,T),V). Hence, the 
particular case of the general theorem of interpolation due to Lions and Magenes [28] 
(see also [35], Chap. Ill, Lemma 1.2), implies that (^r,w) v , = ^^t\ w \ 2 - Taking 
the inner product <-,«;> and using (14) and (16) we conclude that 

^-r-|w| 2 < ^"T-I w l 2 + HHI 2 = - Re < B(w,u),w >< Ci||m|||w| 2 . 

_ Lib _■ (Jib 

Gronwall's inequality now implies 

|w(t)| 2 < |w m | 2 exp (c! J\\u{s)\\ds^ . (34) 

Since w m = u m — v m and u e i 2 ([0, T),V) the equation (32) follows. In particular, 
in the case that u in = v in , it follows that u(t) = v(t) for all t e [0,T]. 

□ 



5. Strong Solutions 

In the analogy to the Navier-Stokes equations we would like to show that the 
shell model equation (9) possesses even more regular solutions under an appropriate 
assumptions. 

Let u m be the solution of the Galerkin system (24) . Taking the scalar product of 
the equation (24a) with Au m , and using Young and Cauchy-Schwartz inequalities 
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we obtain 

~\\u m \\ 2 + v\Au m \ 2 < \(P m B(u m ,u m ),Au m )\ + \(f n ,Au m )\ < 

< Ci||w m || \u m \\Au m \ + i|/| 2 + ^\Au m \ 2 < 

< ^|| u m||2| u m|2 + V _\ Au m \ 2 + 1|/| 2 + j\A U m \ 2 . 

Finally, we get 

By omitting for the moment the term v\Au m \ 2 from the left hand side of the 
equation, we can multiply it by the integrating factor exp(— 2 ^ L J* \u m (s)\ 2 ds) for 
some < to < t to get 

|(||„-||»exp(-^l jf'|«™( fl )| a ifa)) < 2|/|»exp(-^ jfV(,)| a cto) < ^|/| 2 . 

, (35) 

Assuming that the forcing / is essentially bounded in the norm of H for all 
< t < T, namely that / € L°°([0, T], H) and denoting l/^ = ||/|| i ~ ([0 . T]iff) , we 
obtain by integrating equation (35) 

\\u m (t)\\ 2 < e& !> m WV (\\u m (t )\\ 2 + l\ft(t to)) • (36) 

In order to obtain a uniform bound on ||u m (t)|| 2 for all < t < T and for all 
m we need to estimate both the exponent and ||w m (£o)|| 2 - The last quantity is of 
course bounded for to = if we assume that P m u m £ V for all m, however we want 
to put a less strict assumptions. Let us first refine the estimates we made in the 
proof of the existence of weak solutions. Take the scalar product of equation (24a) 
with u m to get 

l^Kf + HKII 2 = (f m , u m ) < ^|/| 2 + ^Kf, 

and using the inequality k\\u m \ 2 < ||u m || 2 we get 

±\u^ + u\\u m f<-^\f\ 2 . (37) 

Integrating from to t, and once again using the fact that / e L°°([0, T], H), we 
get 

v /W(*)ll 2 ^<l« m (0)| 2 + il/&t. (38) 
Jo K \ v 

On the other hand, Gronwall's inequality applied to (37) gives us 

|u m (t)| 2 < \u m {Q)\ 2 e- vk ^ + [ e-^-^^-ds, 

Jo 

and implies 

\u m {t)\ 2 < \u m {0)\ 2 e- k ' vt + (1 - e- H ' vt ), (39) 
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where the last expression yields for every t > 

\u m (t)\ 2 < \u m (0)\ 2 + i^f. (40) 
Using the equation (40) we immediately can evaluate the exponent in (36) by 

e (^/4i« m «r*) < e ^(*-M(i« m (°)i 2 +^i/i 2 ^ (41) 

The last step is to get a uniform bound for the ||« m (io)|| 2 - Let us integrate 
equation (37) from t to t + r. The result, after applying (40) is 

1/ jT^IKfds < \u m (t + r)\ 2 - \u m (t)\ 2 + ^-\f\l < 

<K(0)| 2 + ^|/lL(^; + r 

The standard application of Markov's inequality implies that on every interval of 
length t there exists a time to satisfying 

IK(*o)l| 2 <'0K(O)| 2 + ^l/l^(^; + r)). (42) 

Let us set r < Then, for every t > r and an appropriate to chosen in the 

interval [t — r, t] we get a uniform bound for all m £ N and all t > r, by substituting 
(39) and (42) into (36) 

T\\u™(t)\\* < 2e #^™' 2 +^'^> . ( V(0)| 2 + ^sI/IL) (43) 

To expand the estimate for all times < t < -A^ , we can apply the last inequality 
to the dyadic intervals of the form ^jt+t, ^^r] with r = -^^it+T- Then, if 
* e [fe^7 2' 5 + T ' fef^^l' ^ satisfies r < t < 2r, and for such t inequality (43) implies 

2 ^(i M m (o)i 2 +7 ^i/iD n, ra/ „ xl2 , 3 



t||u™(t)|| a <4e^ vl " Wl T if=*'"~'. ^ K (0)P + ^|/|^j (44) 

Since the last estimate does not depend on fc we can conclude that it holds for all 

< t < 

Finally, passing to the limit in the Galerkin approximation we obtain 

Theorem 4. Let T > and f e L°°([0, T],H). Then if u ln e H, the equation (9) 
possesses a solution u(t) satisfying 

u eL-((o,r],y)nLL((o,r],^(A))fl^([o,T],i/)flL 2 ([o,T],y). 

Moreover, the following bound holds 

sup vklt\\u{t)\\ 2 + sup ||u(t)|| 2 < 

0<t<min{-A-,T} T3~<t<T 



^ a tt^5(I u1 "| 2 + k i 2 \f\1o) (i2\ in\1 , ^ ,,,2 

< 6e fc i" • \ k x \u I + -^2-l/L 



(45) 
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If the initial condition u m £ V , then the solution u(t) satisfies 

u e C([0, T],V)f] L 2 ([0, T],D(A)), (46) 

6. GEVREY CLASS REGULARITY 

6.1. Preliminaries and Classical Navier-Stokes Equation Results. We start 
this section with a definition 

Definition 2. Let us fix positive constants r, p > 0, q > and define the following 
class of sequences 

CO 

(?/ = {«e< a : \e TA "*A«' 2 u\ =^e 2Tfe "fc«K| 2 < oo}. 

n=l 

W^e win say f/iat i/ie sequence u e ^ '« Gevrey regular if it belongs to the class G v T ' q 
for some choice of t, p > and g > 0. 

In order to justify this definition let u — {u\, 112,113, . . . ) and suppose that u n , 
n= 1,2,3,..., are the Fourier coefficients of some scalar function g(x), with corre- 
sponding wavcnumbers k n , defined by relation (2). Then, if u £ G^ q , for some p, q 
and r, then the function g(x) is analytic. The concept of the Gevrey class regularity 
for showing the analyticity of the solutions of the two-dimensional Navier-Stokes 
equations, was first introduced in [20], simplifying earlier proofs. Later this tech- 
nique was extended to the large class of analytic nonlinear parabolic equations in 
[21]. 

For q = we will write G V T — G%'° . Denote by (•, -) PtT the scalar product in G V T 
and 

\ u \p-T = IMIg?, VueGP. 
In the case of q = 1 we will write 

IMIp.t = IM| g p,i, Vu e G p T q , 

and the scalar product in G?' 1 will be denoted by ((•, 

Following the tools introduced in [20] and generalized in [21] we prove the fol- 
lowing: 

Theorem 5. Let us assume that u(0) — u ln e V and the force f £ G? 1 for some 
(7 1 , pi > 0. Then, there exists T, depending only on the initial data, such that for 
every < p < min{pi, log A 1+ 2 V ^ } the following holds 

(i) For every t £ [0, T] the solution u(t) of the equation (9) remains bounded 
m G %)' f or V>W = min{t, cri}. 

(ii) Furthermore, there exists a > such that for every t > T the solution u(t) 
of the equation (9) remains bounded in G? 1 . 

Before proving the theorem, we prove the following estimates for the nonlinear 
term: 

Lemma 6. Let u,v,w £ G V T X for some r > and p < log A ■ Then B(u,v) 
belongs to G V T and satisfies 

\(A 1/2 e rAP/2 B(u,v),A 1/2 e TAP/2 w)\ < C 4 \\u\\ G? ,i ||v|| G? ,i ||w|| G? .i , (47) 
for an appropriate constant C4 > 0. 
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Proof. Let us fix an arbitrary p < log A 1+ 2 V ^ . Observe, that the operator e rAP/2 is 
diagonal in the standard canonical basis {<fi n }%Li 01 H with corresponding eigen- 
values e Tkn , n = 1,2, . . . , where k n = A™ for some A < 1+ 2 V ^ . Therefore, 

\{A 1 ' 2 e rA, "*B{< U ,v),A 1 / 2 e rA " 2 w)\ < 



< 



n=l 
1,2 



ak n+1 k 2 n v n+2 u* n+1 w* n + bklv n+1 w* n u* n _ 1 



+ ak n k n - 1 w*u n - 1 v n - 2 + bk*k n -iw*v n -iu n - 2 



< 



< sup|e rfc "fc„w„| e Tk " 

n>l f"~i 



+ 



ak n k n -iU n -iV n - 2 



ak n+ ik n v n+2 u* n+1 
bk n k n -iu n - 2 v n -i 



+ 



+ 



< C 4 ||u;|| G¥ ,i||u|| G? ,i||w|| G? ,i, 

(48) 

where we can choose C4 = (|a|(A~ 2 + A 2 ) + |6|(1 + A 2 )). The last inequality is 
the result of the Cauchy-Schwartz inequality and the fact, that e T ( A "- A " ) ; 
e 7-(A n -A n+1 -A n 1 ) an( j g r(A"-A" X -A" 2 ) are j egg than or e q ua i to 1 for all t > and 

our specific choice of A. 

□ 

Now we are ready to prove the Theorem: 

Proof. (Of Theorem 5) Let us take an arbitrary < p < min{p!,log A 1+ 2 } an d 
denote ip(t) = min(i, <j\). For a fixed time t take the scalar product of the equation 
(9a) with Au(t) in G^( t ) to get 

(u\t),Au(t)) pMt) +v\Au(t)\l Mt) = 

= (f,Au(t)) pMt) - (B(u(t),u(t)),Au(t)) pMt) . (49) 
The left hand side of the equation could be transformed in the following way: 
{u'{t),Au{t)) P Mt) = (e^ AP/2 u'(t),e^ AP/2 Au(t)) = 

= (A 1 / 2 (^(')^ /2 u (t))' - <S(t)A^ 1 V 2 e*W A ''\(t),AV 2 e*W AP '\(t» = 



1 d 



¥>(*) 



> 



2 dV 
ld_ 

2dt' 

> lj t \Ht)\\ 2 pMt) - %Au(t)\l v(t) l\\u(t)\\ 2 pMt) . 

The right hand side of the equation (49) can be bounded from above in the usual 
way using Lemma 6, Cauchy-Schwartz and Young inequalities 

\(f,Au(t)) pMt) - {B(u{t),u{t)),Au{t)) pMt) \ < 



> 9^ll«(*)llp lV (t) ~ I^WU( t )IKi)IU( t ) > 



(50) 



< 



\f\p, v (t)\Au{t)\ p ^ [t) + C 4 ||u||p )¥ , (t) 



< 



<-|/lLw + 4l^(*)lLw+^ll^llLw 



(51) 
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Combining (50) and (51), the equation (49) takes the form 

j t \Ht)\\ 2 pMt) + »\Au(t)\ 2 pMt) < 2 -\f\i Mt) + Inmirtt) + c4< M t) * 

<l\f\l v{ t) + C 5 + C 4 \\u\\l Mt) . 

Denoting y(t) = 1 + \\u(t)\\ 2 p , v(t ) and K = §l/lp, v (t) + <3> + C 4 2/3 we get an 
inequality 

It implies 

y(t) = l + |Kt)||^ (t) <2 tf (0) = 2(l + ||« in || 2 ), 

for 

t < T\ (o"i , ||u*1 , |/U (t) ) = 2 -^y- 1/2 (Q) = 2 -^(i + W^W*)- 1 '*. 

According to Theorem 4, the solution u(t) of the equation (9) remains bounded 
in V for all time if we start from u m G V. Moreover, from the relation (46) we 
conclude that there exists a constant M, such that for all t > 

||u(t)|| < M. 

Hence we can repeat the above arguments starting from any time t > and find 
that 

\e^ AP/2 A 1 / 2 u(t)\<2 + 2M 2 , (52) 
for all t>T 2 = 2^2(1 + M 2 )' 1 / 2 , and the Theorem holds. 

□ 

6.2. A Stronger Result. We prove the following stronger version of Theorem 5; 
it is not known whether this holds for the Navier-Stokes equations. 

Theorem 7. Let us assume that u(0) = u m € H and the force f € G?\ for some 
o~\,P\ > 0. Then, there exists T, depending only on the initial data, such that for 
every < p < min{pi, log A 1+ 2 V ^ } the following holds 

(i) For every t G [0, T] the solution u{t) of the equation (9) remains bounded 
in G ^(t)> / or V'(i) = min{i,eri}. 

(ii) For every t > T there exists a > such that the solution u{t) of the 
equation (9) remains bounded in G v a . 

Before proving the theorem, we obtain the following estimates for the nonlinear 
term: 

Lemma 8. Let u, v G G V T , w G G? A for some t > and p < log A 1+ 2 V ^ ■ ■ Then 
B(u,v) belongs to G V T and satisfies 

\{e TAP/2 B(u,v),e TAP/2 w)\ < C 6 ||u|| G? ||w|| G? ||u;|| G? ,i, (53) 
for an appropriate constant C§ > 0. 
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Proof. Let us fix an arbitrary p < log A 1+ 2 V ^ . Then we note that the operator 

e rAP/2 is diagonal with eigenvalues e T ^" , n = 1, 2, . . . , where fc„ = A" for A < 1+ 2 V ^ ■ 
Hence, 



(e^ /J 5(u,t)) 1 e TX ' /J «))| < 



OO 
71=1 



a/j„ + i?i n+2 w* +1 w* + 



ck n — 



< sup|e Tfc,, fc„w„| ^ e T,! 



n>l 



a\u n+2 v* n+1 



< 



+ 



CA 1 V n -lU n - 2 



< 



< C 6 ||ti;|| G? ,i||«|| G? ||«|| G P, (54) 

where we can choose Cq = (Aa + 6 + A _1 c). The last inequality is the result of the 
Cauchy-Schwartz inequality and the fact, that e T ( A ~ A + ~ A + \ e T ( A ~ A + ~ A ) 
and e T ( A "~ A " ~ A " ) are less than or equal to 1 for all r and the specific choice of 
A. 

□ 

Now we are ready to prove the Theorem: 

Proof. (Of Theorem 7) Let us take an arbitrary < p < min{pi,log A 1+ 2 } and 
denote <p(t) = min(i, a{). For a fixed time t take the scalar product of the equation 
(9a) with u(t) in GP(t) to get 

(u'(t),u(t)) pMt) + i/||«(t)||p iV , (t) = 

= (/,«(*) W) - (B(u(t),u(t)),u(t)) pMt) . (55) 
The left hand side of the equation could be transformed in the following way: 
(u'(t),Au{t)) p , v(t) = (e^ AP/2 u'(t),e^ AP/2 u(t)) = 

= ((e^ AF/2 u(t)y - p'(t)Av/ 2 e^ AP/2 u(t),e^ AP/2 u(t)) = 

^l«(*)lp, v (t)-^W(^ /2 «W.«WW(t)> 



2 dt 
1 d 



2 dt 
1 d 



> 



(56) 



>^l«(*)l^(t)-ill«(*)llp, V (t)-JK*)lp, V ( t )- 

The right hand side of the equation (55) could be bounded from above in the 
usual way using Lemma 8, Cauchy-Schwartz and Young inequalities 



!(/>«(*) W(t) ~ (B(u(t),u(t)),u(t)) pMt 



< 



< \f\p, v (t)\u(t)\ PlV (t) +C6\\u\\ p , ip(t) \u\l Mt) 



z ii/Cw + si«(*)i^w + iNiU + (57) 



< 
r 2 



18 



P. CONSTANTIN, B. LEVANT, AND E. S. TITI 



Combining (56) and (57), the equation (55) takes the form 

^Ht)\lrtt) + v \\Mt)\\l, v <t) - 

^l/&) + ^K*)lL( t ) + (^ + ^)l<,( t )< 

<l\K% + C 7 \u\$, v( t), 

for some positive constant C-j. 

Denoting y(t) = 1 + \u{t)\ 2 pv{t) and K = §|/|p^ (t) + C 7 we get an inequality 

V < Ky 2 , 

It implies 

y{t) = l + \\u{t)\\ 2 p ^ t) <2y{Q)^2{l + \u m \ 2 ) 1 

for 

t < ^(tn, |« <n |, |/U (t) ) = ^y-\o) = 2^(1 + lu*"! 2 )" 1 . 

According to Theorem 2 the solution of the equation (9) remains bounded in H 
for all time, if we start from u m £ H. Moreover, from the estimate (23), it follows 
that there exists a constant M, such that for all f > 

|u(t)| < M. 

Hence we can repeat the above arguments starting from any time t > and find 
that 

AV 2 u(t)\ < 2 + 2M 2 , (58) 
for alii > T 2 = ^(1 + M 2 ) -1 , and the Theorem holds. 

□ 

7. Global Attractors and Their Dimensions 

The first mathematical concept which we use to establish the finite dimensional 
long-term behavior of the viscous sabra shell model is the global attractor. The 
global attractor, A, is the maximal bounded invariant subset of the space H. It 
encompasses all of the possible permanent regimes of the dynamics of the shell 
model. It is also a compact subset of the space H which attracts all the trajectories 
of the system. Establishment of finite Hausdorff and fractal dimensionality of the 
global attractor implies the possible parameterization of the permanent regimes of 
the dynamics in terms of a finite number of parameters. For the definition and 
further discussion of the concept of the global attractor see, e.g., [14] and [36]. 

In analogy with Kolmogorov's mean rate of dissipation of energy in turbulent 
flow we define 

the mean rate of dissipation of energy in the shell model system. (•) represents 
the ensemble average or a long-time average. Since we do not know whether such 
long-time averages converge for trajectories we will replace the above definition of 
e by 

1 /"* 

e = v sup limsup- / 2 <is. 

u in eA t^oo t Jo 
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We will also define the viscous dissipation length scale Id- According to Kol- 
mogorov's theory, it should only depend on the viscosity v and the mean rate 
of energy dissipation e. Hence, pure dimensional arguments lead to the definition 

1/4 

Id 



which represents the largest spatial scale at which the viscosity term begins to 
dominate over the nonlinear inertial term of the shell model equation. In analogy 
with the conventional theory of turbulence this is also the smallest scale that one 
needs to resolve in order to get the full resolution for turbulent flow associated with 
the shell model system. 

We would like to obtain an estimate of the fractal dimension of the global attrac- 
tor for the system (9) in terms of another non-dimensional quantity - the generalized 
Grashoff number. Suppose the forcing term satisfies / G L°°([0, T],H) and denote 
11/11 i°°([o oo) H) ~ Then we define the generalized Grashoff number for our 

system to be 

The generalized Grashoff number was first introduced in the context of the study of 
the finite dimensionality of long-term behavior of turbulent flow in [15]. To check 
that it is indeed non-dimensional we note that ||/|| £ ooq t , h ^ has the dimension of 

^2 , where L is a length scale and T is a time scale. k\ has the dimension of j- and 
the kinematic viscosity v has the dimension of 4|r. In order to obtain an estimate 
of the generalized Grashoff number, which will be used in the proof of the main 
result of this section, we can apply inequality (38) to get 



lim 

t 



nsupi f \\u(s)\\ 2 ds< !& = v 2 k\G 2 . (60) 
— >oo t Jq v k^ 

Theorem 9. The Hausdorff and fractal dimensions of the global attractor of the 
system of equations (9), dff(A) and dp (A) respectively, satisfy 

d H (A) < d F (A) < log A (f-) + ilog A fc.iX 2 - 1)). (61) 



In terms of the Grashoff number G the upper bound takes the form 
d H (A) < d F (A) < log A G 1 ' 2 + l - log A (c^X 2 - 1)) . 



(62) 



Proof. We follow [6] (see also [4], [5], [36]) and linearize the shell model system of 
equations (9) about the trajectory u(t) in the global attractor. In the Appendix B 
we show that the solution u{t) is diffcrentiable with respect to the initial data, 
hence the resulting linear equation takes the form 

^ + vAU + B (t ) U = (63a) 
U(0) = U m , (63b) 

where B (t)U = B(u(t),U(t)) + B(U(t) 7 u(t)) is a linear operator. In order to 
simplify the notation, we will denote 

A(t) = -vA - B (t). 
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Let Uj(t) be solutions of the above system satisfying Uj(0) — Uf 1 , for j = 
1,2, ...,N. Assume now that U\ n , U™, ■ ■ ■ , U™ are linearly independent in H, 
and consider QN(t) - an JT-orthogonal projection onto the span {Uj{t)}^ =l . Let 
{(Pj(t)}^ =1 be the orthonormal basis of the span of {Uj{t)}^ =1 . Notice that ipj € 
D(A) since spaniU^t), U N {t)} C D(A). 

Using the definition of A(t) and the inequalities of the Proposition 1 we get 

N 

Re(Trace[A(t) o Q N (t)}) = Re ^(A(i)v?j , <pj) = 

3=1 

N 

= E —^11 fi\\ 2 - MBfaMt)), fj) < 

N 

< -vY^itf + \{B{<p i ,u(t)),y j )\) < 

3 = 1 

\2JV _ i 

< - 1 yk 2 \ 2 ^- Y + C\N\\u(t)\\. (64) 

According to the recent results of [4] (see also [5], [6], [36]), if N is large enough, 
such that 

1 /"' 

limsup- / Re(Trace[K(t) o Q N (s)])ds < 0, 

t— tOO t Jq 

then the fractal dimension of the global attractor is bounded by N. We need 
therefore to estimate ./V in terms of the energy dissipation rate. Using (64), the 
definition of e we find that it is sufficient to require N to be large enough such that 
it satisfies 

\ 2N - 1 / 1 /•* \ 1/2 / f\ 1/2 

^o 2 A 2 ^^ r >C 1 iV^ limsup- \\u(t)\\ 2 j =CiN\j) . (65) 

Finally, tv~ 3 = l^ 4 implies 

which proves (61). Applying the estimate (60) to the inequality (64) we get the 
bound (62) in terms of the generalized Grashoff number. 

□ 



8. Determining Modes 

Let us consider two solutions of the the shell model equations u, v corresponding 
to the forces /, g e L 2 ([0, oo), H) 

du 

— + vAu + B(u, u) = f, (66) 
dv 

— +vAv + B(v,v)=g. (67) 

We give a slightly more general definition of a notion of the determining modes, 
than the one that is introduced previously in literature [16] (see also [14], [15], [23] 
and references therein). 
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Definition 3. We define a set of determining modes as a finite set of indices 
McN, such that whenever the forces f, g satisfy 

\f(t)-g{t)\^0, ast^oo, (68) 

and 

\ u n(t) - v n (t)\ 2 -> 0, ast^oo (69) 

it follows that 

\u(t) -v(t)\ -> 0, ast^oo. (70) 

TTie number of determining modes N of the equation is the size of the smallest such 
a set M. 

We would like to recall the following generalization of the classical Gronwall's 
lemma which was proved in [23], [24] (see also, [14]). 

Lemma 10. Let a = a(t) and (3 = (3{f) be locally integrable real-valued functions 
on [0, oo) that satisfy the following condition for some T > 0: 



1 f t+T 

- / a(r)dr > 0, (71) 
1 f t+T 

limsup— / a (r)dr < oo, (72) 
t^oo T J t 

i r t+T 

lim - / f3+(r)dT = 0, (73) 

where or (t) = max{— a(t), 0} and f3 + (t) = max{/3(i), 0}. Suppose that £ — £(i) is 
an absolutely continuous nonnegative function on [0, oo) that satisfies the following 
inequality almost everywhere on [0, oo): 

§ + «f</3. (74) 

Then £(t) —^0 as t — > oo. 

A weaker version of the above statement was introduced in [15]. In fact this 
weaker version would be sufficient for our purposes. 

Theorem 11. Let G be a Grashoff number defined in (59). Then, 

(i) The first N modes are determining modes for the shell model equation 
provided 

N> l -\og x {C 1 G). (75) 

(ii) Let u,v be two solutions of the equations (66) and (67) correspondingly. 
Let the forces f,g satisfy (68) and integer N be defined as in (75). If 



lim inf 

t— >oo T 



lim |ujv(t) — VN(t)\, and lim \uN-i(t) — VN-i(t)\ = 0, 



then 



lim \Q N u(t) - Q N v(t)\ =0. 

i— >-oo 
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Proof. Denote w = u — v, which satisfies the equation 

^ + vAw + B(u, w) + B{w,v) = f - g. (76) 

Let us fix an integer m > 1 and define P m : H — > H to be an orthogonal projection 
onto the first m coordinates, namely, onto the subspace of H spanned by m vectors 
of the standard basis {4>i\YLi- Also denote Q m = I — P m . 

It is not hard to see that the second part of the Theorem implies the first state- 
ment. Therefore it is enough for us to prove that if 

lim |iojv(t)|, and lim |iojv_i(t)| = 0, 



then for every m > | log A (CiG) 



lim |Q m tu| = 0. 



Multiplying the equation (76) by Q m w in the scalar product of H we get 

~\Q m w\ 2 + v\\Q m w\\ 2 = 

= Re(f - g, Q m w) - Re(B(u, w), Q m w) - Re(B(w, v),Q m w). (77) 

For the left hand side of the equation we can use the inequality k 2 n+1 \Q m w\ < 
\\Q m w\\ . For the right hand side of the equation we observe that 

oo , 

Re(B(u, w), Q m w) = Im ^ ( ak n+1 w n+2 u* n+1 w* n + bk n w n+1 w* n u* n _ 1 + 

n— m+1 ^ 

+ ak n ^ 1 W* n U n ^lW n ^ 2 + i)fc„-lW> n -lt(„-2 I < 



< 



ak m W m ^_iU m W m ^i + ak m J r iW m ^2 u m+l w rn + bk m W m ^_iW m U m —i 



Moreover, 



OO / 

Re(B(w, v),Q m w) =im )^ ( ak n+ iV n+ 2W* n+1 w* n + bk n v n+1 w* n w* n _ v 

n — ra -I- 1 ^ 



n—rn-\-l 

+ ak m+1 w* m+2 w m+1 v m + 



+ bk m w* m+1 v m w m -i 



+ Ci\\v\\ \Q m w\ 2 . 



Finally denoting £ = |Q m u>| 2 we can rewrite the equation in the form 



ati < a, 



where 



2 dt 

a{t) = uk 2 n+1 -C 1 \\v{t)\\, 
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and 

/?(*) = (f - 9, QmW) + 



+ 



ak m w' m+1 u m w m -i + ak m+1 w* n+2 u m+1 w m + bk m w 
ak m w* m+1 w m v m - 1 + ak m+1 w m+2 w m+1 v m + bk m w* m+1 v m w m -i 



To see that (3(t) satisfies condition (73) of Lemma 10 it is enough to show that 
(3(t) — ► as t — > oo. It is clearly true, if we assume that the forces /, g satisfy (68) 
and |w m (t)|, |u; m _i(i)| — > as t — > oo. Moreover, it is easy to see that a(t) satisfies 
the condition (72). Hence, in order to apply Lemma 10 we need to show that a(t) 
also satisfies the conditions (71), namely, that 

ft+T / ! r t+T \ 1/2 

\v(t)\\ 2 (1t) >0. (78) 



i r +T ( i r 

liminf - / a(r)dT > uk 2 m+1 - C x limsup - / 



To estimate the last quantity, we can use the bound (60) to conclude that if m 
satisfies 

A 2m > dG, 

then 

|Q m to| — » as t — > oo. 

□ 

Remark 1 . The existence of the determining modes for the Navier-Stokes equations 
both in two and three dimensions is known (see, e.g., [7] [14], [24] and reference 
therein). However, there exists a gap between the upper bounds for the lowest 
number of determining modes and the dimension of the global attractor for the 
two-dimensional Navier-Stokes equations both for the no-slip and periodic boundary 
conditions. Our upper bounds for the dimension of the global attractor and for the 
number of determining modes for the sabra shell model equation coincides. Recently 
it was shown in [22] that a similar result is also true for the damped-driven NSE 
and the Stommel-Charney barotropic model of ocean circulation. 

9. Existence of Inertial Manifolds 

In this section we prove the existence of a finite-dimensional inertial manifold 
(IM). The concept of inertial manifold for nonlinear evolution equations was first 
introduced in [17] (see also, e.g., [5], [8], [9], [18], [34], [36]). An inertial manifold is a 
finite dimensional Lipschitz, globally invariant manifold which attracts all bounded 
sets in the phase space at an exponential rate and, consequently, contains the global 
attractor. In fact, one can show that the IM is smoother, in particular C 1 (see, e.g., 
[10], [32], [34]). The smoothness and invariance under the reduced dynamics of the 
IM implies that a finite system of ordinary differential equations is equivalent to the 
original infinite system. This is the ultimate and best notion of system reduction 
that one could hope for. In other words, an IM is an exact rule for parameterizing 
the large modes (infinite many of them) in terms of the low ones (finitely many of 
them). 

In this section we use Theorem 3.1 of [18] to show the existence of inertial 
manifolds for the system (9) and to estimate its dimension. Let us state Theorem 
3.1 of [18] in the following way (see also [36], Chapter VIII, Theorem 3.1) 
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Theorem 12. Let the nonlinear term of the equation (81) R(-) be a differentiable 
map from D{A) into D{A 1 ^ 13 ) satisfying 

\R'(u)v\<C 8 \Au\\A^v\ (79a) 

lA^PR'^vl < C 9 \Au\\Av\, (79b) 

for some < [3 < 1, for all u, v G D(A) and appropriate constants C 8 ,C g > 0, 
which depend on physical parameters a, b, c, v and f . 

Let b > and I <G [0, 1) be two fixed numbers. Assume that there exists an N, 
large enough, such that the eigenvalues of A satisfy 

where 

K 2 = ^\A 1 - f3 f\ + 26C 9P , 

p is the radius of the absorbing ball in D{A), whose existence is provided by Propo- 
sition 14- 

Then the equation (81) possesses an inertial manifold of dimension N which is 
a graph of a function $ : PnH — > QnD{A) with 

b= sup 
{ P eP N H} 

and 

l = \A^( Pl ) - ^( P2 ))\ ^ 

{pi,P2ePNH, Pl7 t P2 } \A(Pi-P2)\ 
is the Lipschitz constant o/$. 

In order to apply the Theorem, we need to show that the sabra model equation 
satisfies the properties of the abstract framework of [18]. Let us rewrite our system 
(9) in the form 

du 

— + vAu + R(u) = 0, (81) 

where the nonlinear term R(u) = B(u,u) — f. 

First of all, we need to show that the nonlinear term R(u) of the equation (81) 
satisfies (79) for [3 = 0. Indeed, and based on (15) R(u) is a differentiable map 
from D(A) to Z^A 1 / 2 ) = V, where 

R'(u)w = B(u,w) +B(w,u), Vu,weD(A). 

Moreover, the estimates (79) are satisfied according to Proposition 1, namely 

\R'(u)v\ = \B{u,v) + B{v,u)\ < (d+C 2 )HM < (Cl ± ° 2) \Au\\v\, 

ki 

where the last inequality results from the fact that fci < \Au\ and the constant 

C 8 = £l+£A = (| a |(A- 1 +A) + |h|(A^ + l)) + (2|a|+ 2 A| b |) ^ ^ ^tion, Qnc ^ proye ag ^ 

Proposition 1 that 

\AR'(u)v\ < \AB(u,v)\ + \AB(v,u)\ < C 9 \\u\\\Av\ < C 9 \Au\\Av\, 
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where here again the last inequality results from the fact that < \Au\ and 

the constant Cg = Cg/fci- 

Finally, because of the form of the wavenumbers k n (2), the spectral gap condi- 
tion (80) is satisfied for (3 = 0. Therefore, we can apply the Theorem 12 to equation 
(9), and conclude that the sabra shell model possesses an inertial manifold. 

9.1. Dimension of Inertial Manifolds. Let us calculate the dimension of the 
inertial manifold of the the sabra shell model equation for the specific choice of 
parameters 

= 0, 1 = 1, b = p, 

where p is the radius of the absorbing ball in the norm of the space D (A) (see 
Proposition 14). In that case, the conditions (80) take the form 

k N+1 -k N >max^K 2 ,^-j, (82) 

and 

Corollary 13. The equation (9) possesses an inertial manifold of dimension 

AT > max {log^,log AII ^}. (83) 

10. Conclusions 

We have established the global regularity of the sabra shell model of turbu- 
lence. We shown analytically that the shell model enjoys some of the commonly 
observed features of real world turbulent flows. Specifically, we have established 
using the Gevrey regularity technique the existence of an exponentially decaying 
dissipation range, which is consistent with the observations established in [33] for 
the GOY model using asymptotic methods. Moreover, we have provided explicit 
upper bounds, in terms of the given parameters of the shell model, for the number 
asymptotic degrees of freedom. Namely, we presented explicit estimates for the 
dimension of the global attractor and for the number of determining modes. In the 
forthcoming paper we will investigate the lower bounds on the Hausdorff dimension 
of the global attractor. 

Finally, we have shown that the shell model possess a finite dimensional inertial 
manifold. That is, there exists an exact rule which parameterizes the small scales 
as a function of the large scales. The existence of inertial manifolds is not known 
for the Navier-Stokes equations. It is worth mentioning that the tools presented 
here can be equally applied to other shell models. 
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Appendix A. Dissipativity in D{A) 

In this section we prove that the equation (9) is dissipative in different norms 
and that its solutions are real analytic with respect to the time variable with values 
in D(A). 

First, we start with the space if. Consider the relation (39), and by letting 
t^oowe obtain 

limsupKi)! < = Gvki, (84) 

where |/|oo = |j/|| L ^^ Q ^ H \ and G is the generalized generalized Grashoff number, 
defined in (59). 

The existence of the absorbing ball in the space V readily follows from Theorem 4. 
However, the exact formula for the radius of the absorbing ball is much more 
involved. 

We already mentioned that if we assume that the forcing / belongs to some 
Gevrey class for all t > 0, then the existence of absorbing balls for solutions of the 
equation (9) in the norms of V and D(A) follows from the fact that in that case 
all the functions belonging to the global attractor arc in some Gevrey class and 
hence are bounded in all weaker norms. In the current section we will show that 
the strong assumption on the force / to be in the Gevrey class could be dropped. 
In particular, we will assume that / G if is time independent. Our proof follows 
[19] (see also [5], ch. 12). 

Fix m G N, let P m be the projection of if onto the first m coordinates and 
denote 

P m u = u m . 

First, let us complexify our equation and all the relevant spaces and operators. Re- 
call that the Galerkin approximating system (24) of order m for equation (9) is an m 
- dimensional system of ordinary differential equations, with analytic nonlinearity, 
and hence its solutions are locally complex analytic functions. 

Let us consider the complex time variable t = se l9 , 9 G (—§,§) and s G R+. 
Then 

^-\\u m (se ie )\\ 2 = ^-{u m (se ie ),Au m {se 10 )) = 2Re(e t9 ^-u m 1 Au m ). 

do (Xb (XL 

Then, by taking the scalar product of (24a) with Au m , multiplying by e ld and 
taking the real part we obtain 

~\\u m (se ie )\\ 2 + v cos 9\Au m (se ie )\ 2 = 

_ (IS 

= Re (e l9 {f : Au rn ) - e w {P m B(u m , u m ), Au m )j < 

< ^\Au™{se*)\ 2 + -g- + -^-\\u-{ S e^)\\\ 
2 v cos a v cos a 

where the last line follows from (18) and Young's inequality. Finally we deduce 

^-\\u m ( S e* 9 )\\ 2 + uco S 0\Au m (se t9 )\ 2 < + ^-\\u m {se M )\\\ 

as v cos v cos 

Now we are able to derive the bound 

|| U m (t)|| 2 < 2(|| u m (0)|| 2 + 1) < 2(|K"f + 1), (85) 
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provided t — se l ° satisfies 



*< 2^(||« in || 2 + l)-\ (86) 



where K = -^0^ + -^L. 

v cos v cos 

The bounds (85) and (86) show that for every m£N, the functions u m (t) : C — > 
C m are analytic in the domain 

D = D{v,\n\f\) = {t = se ie : |0| < J, < a < ^(IkT + l)" 1 }. 

Let 7 be a small circle contained in the domain D. Then according to Cauchy 
formula 

d 



If u m (z) , 
Using the relation (85) we get 



dt 



\\ju m \\<^{\\u™{0)\\ 2 + l) 112 , 

where r 7 is the radius of 7. Let M C D be a compact subset of _D and denote 
r M = dist(M,dD), then 

||* m ||<H^(|| u ™ ( o)f + 1)V2 ; (87) 

for all t e M and all m e N. 

Finally comparing the H norms of both sides of the equation (24a), using the 
relations (12), (87) and the triangle inequality and we get 

u\Au m {t)\<\ju m \ + \f{t)\ + \B{u m ,u m )\< 

<k- 1/2 \\ju m \\+\f(t)\ + \B{u m ,u m )\< 

< ^L(i| U ™(o)n 2 + 1) 1 / 2 + \f(t)\ + cAvrm^m, 

for all t £ M. Since all the summands at the right hand side are bounded (see the 
proofs of Theorem 2 and Theorem 4) and those bounds do not depend on m we 
can pass to the limit in m, using Vitali's convergence theorm for complex analytic 
functions, concluding that there exists a constant E = E[y, |/|, ||u m || , A, fc , a, b, c), 
independent of m, such that 

\Au(t)\ < E, 

for all t £ M. Moreover, for all t £ M the solution u{t) is a complex analytic 
function with values in D(A). Repeating this procedure starting instead of from 
some tx E M we conclude that is uniformly bounded for all t > 0. 

Proposition 14. Let f £ H and u m £ V . Then the solution of the equation (9) 
is an analytic function, with respect to the time variable with values in D(A), and 
that it possesses absorbing balls in the spaces H,V and D(A). 
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Appendix B. Differentiability of the semigroup with respect to 

INITIAL CONDITIONS 

Theorems 2 and 3 shows in fact, that the initial-value problem (9) is well posed. 
This allows us to define the semigroup S(t), i.e., the one parameter semigroup 
family of operators 

S(t) :u m eH^ u(t) e H, 
which are bounded for almost all t > 0. According to the Theorem 3, the mapping 
S(t) : H —* H is Lipschitz continuous with respect to the intial data. The purpose 
of this section is to show that this mapping is also Frechet differentiable with respect 
to the initial data (see, e.g., [36] (ch. VI, 8)). 

First of all we wish to linearize the non-linear term. Let us fix T > and let 
u, v e L°°([0, T],H)n L 2 ([0, T], V) be solutions of the equation (9) with the initial 
conditions u(0) = u m and v(0) = v m . Then for almost every t G [0, T] we can write 

B(u(t), u(t)) - B(v(t),v(t)) = B (t)(u(t) - v(t)) + B^t, u - v), 

where Bq (t) : H — > V is a linear operator defined as 

B (t)w = B(u(t),w(t)) + B(w(t),u(t)), 

and 

Si (*,«;) = -B(w(t),w(t)). 
Let us consider a solution U (t) of the linearized equation 

^ + vAU + B (t ) U = (88a) 
U{0) = U in =u in -v in , (88b) 
satisfying Ue L°°([0,T],H)n L 2 ([0,T],V). Denote 

f = u — v — U. 

Using the fact that u, v satisfy equation (9) and U is the solution of (88) we can 
directly check that ip satisfies 

^ + vAip + B (t)if = -Bi(t, u(t) - w(t)), (89) 
and <p(0) = 0. Multiplying by if we get 

\j t \^\ 2 + "M 2 + (Bo(t)<p,<p) - - (Bi(t,u(t) - w(t)), tp) . (90) 
According to the Proposition 1 we can rewrite the last equation as the inequality 

\j t W\ 2 + HMI 2 < c'\p\ ■ \u\ ■ \\f\\ +c 9 \u- v\ 2 ■ \\f\\ . 

The fact that u E L°°([0, T], H) together with Young's inequality imply 

Our next step is to use Gronwall's inequality to get 

|^)|2< 9l f ' e^( T s)\u{ s )-v{stds, 
v Jo 

for every t G (0,T]. Once again applying Theorem 3 we finally get 

W)\ 2 <C w \u m -v ln \\ 
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2„4K 



where < C\q = — ^ — J e » ( T s ^ds and K is the constant from the statement 
of the Theorem 3. It follows that 

Ht)-vit)-u{t)\ < Q< 

\u m - v m \ ~ ~ 

as \u m — v m \ tends to 0. That exactly means that U(t) is a differential of S(t) with 
respect to u m e H. 
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